N94-21477 


TWO-LEVEL SCHWARZ METHODS FOR 
NONCONFORMING FINITE ELEMENTS 
AND DISCONTINUOUS COEFFICIENTS’ 


Marcus Sarkis 
Courant Institute 
New York, NY 


f. 


SUMMARY 


Two-level domain decomposition methods are developed for a simple nonconforming 
approximation of second order elliptic problems. A bound is established for the condition number 
of these iterative methods, which grows only logarithmically with the number of degrees of freedom 
in each subregion. This bound holds for two and three dimensions and is independent of jumps in 
the value of the coefficients. 


INTRODUCTION 

The purpose of this paper is to develop domain decomposition methods for second order elliptic 
partial differential equations approximated by a simple nonconforming finite element method, the 
nonconforming P\ elements. We consider a variant of a two-level additive Schwarz method 
introduced in 1987 by Dryja and Widlund [1] for a conforming case. In these methods, a 
preconditioner is constructed from the restriction of the given elliptic problem to overlapping 
subregions into which the given region has been decomposed. In addition, in order to enhance the 
convergence rate, the preconditioner includes a coarse mesh component of relatively modest 
dimension. The construction of this component is the most interesting part of the work. Here we 
have been able to draw on earlier multilevel studies, cf. Brenner [2], Oswald [3], as well as on recent 
work by Dryja, Smith, and Widlund [4]. Our main result shows that the condition number of our 
iterative methods is bounded by C (1 + log(H/h), where H and h are the mesh sizes of the global 
and local problems, respectively. We also note that this bound is independent of the variations of 
the coefficients across the subregion interfaces. 

The face based and the Neumann-Neumann coarse spaces, that we are introducing, have the 
following characteristics. The nodal values are constant on each edge (or face) of the subregions 
and the values at the other nodes are given by a simple but nonstandard interpolation formula. 
Thus the value at any node in the interior of a subregion is a convex combination of three (or four) 
values given on the boundary, in case of triangular (or tetrahedral) substructures. We note that an 
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important difference between nonconforming and the conforming case is that there are no nodes at 
the vertices (or wire basket) of the subregions. 

We note that ideas similar to ours have been used recently in other studies of domain 
decomposition methods for nonconforming elements; cf. Cowsar [5,6] and Cowsar, Mandel and 
Wheeler [7], In particular, an isomorphism similar to ours was independently introduced by 
Cowsar. We point out that by using these isomorphisms, we can analyze any nonconforming 
version of domain decomposition methods which have already been analyzed for conforming cases. 
In this paper, we focus on the case where there are great variations in the coefficients across 
subdomains boundaries for both two and three dimensions. We define and analyze new coarse 
spaces and obtain condition numbers with just one log factor. 

A short version of this paper was entered into Copper Mountain student competition in 
mid-December 1992. The present paper is a slight modification of a technical report [8]. 

DIFFERENTIAL AND FINITE ELEMENT MODEL PROBLEMS 


To simplify the presentation, we assume that Q is an open, bounded, polygonal region of 
diameter 1 in the plane, with boundary d£l. In a separate section, we extend all our results to the 
three dimensional case. 

We introduce a partition of f2 as follows. In a first step, we divide the region Q into 
nonoverlapping triangular substructures *, i = 1, • • • , N. Adopting common assumptions in finite 
element theory, cf. Ciarlet [9], all substructures are assumed to be shape regular, quasi uniform and to 
have no dead points; i.e. each interior edge is the intersection of the boundaries of two triangular 
regions. We can show that the theory also holds if we choose nontriangular substructures, where 
the boundary of each substructure is a composition of several curved edges, and each curved edge is 
the intersection of two substructures. Naturally, we need assumptions related to the quasi 
uniformity and nondegeneracy of this partition. Initially, we restrict our exposition to the case of 
triangular substructures since the main ideas are seen in this case. This partition induces a coarse 
mesh and we introduce a mesh parameter H := max{ifi, • • ■ , H^} where Hj is the diameter of flj. 

We denote this triangulation by T H . Later, we extend the results to nontriangular substructures. 

In a second step, we obtain the elements by subdividing the substructures into triangles in such 
a way that they are shape regular, and quasi uniform. We define a mesh parameter h as the 
diameter of the smallest element and denote this triangulation by T h . Similarly, we assume the 
triangulation T h does not have any dead points. 

We study the following selfadjoint second order elliptic problem: 

Find u e Hq(Q), such that 

a(u,v) = f(v), V v € H^(Q) , (1) 


544 


where 


a(u, v) = J a(x) Vtt • Vv dx and f(v ) = J^fvdx for / G L 2 . 

We assume that g(x) > a > 0 and that it is a piecewise constant function with jumps occurring 
only across the substructure boundaries. This includes cases where there is a great variation in the 
value of the coefficient a(x). We remark that there is no difficulty in extending the analysis and the 
results to the case where a(x) does not vary greatly inside each substructure. 

Definition 1 The nonconforming P\ element spaces (cf. Crouzeix and Raviart. [10]) on the h-mesh, 
and H-mesh is given by 

V h := {u | v linear in each triangle T G T h , 

v continuous at the midpoints of the edges of T h , and 
v = 0 at the midpoints of edges ofT h that belong to 

and 

V H := {u | v linear in each triangle T G T H , 

v continuous at the midpoints of the edges ofT H , and 
v = 0 at the midpoints of edges of T H that belong to cK)} . 

These spaces are nonconforming; in fact V H V h and V h <£. Hq(CI). 

Let S be a region contained in Cl such that does not cut through any element. Denote by 
and T h |£ the space V h and the triangulation T h restricted to S, respectively. 

Given u G Vjj, we define the discrete weighted energy semi norm by: 

MV (S) :=4(u,u). ( 2 ) 

a,h x 1 

where . 

a^(u,v) = ^2 / a(x)Vu-Vudx. (3) 

TeT» lt jT 

In a similar fashion, we define the inner product Oq(u,v) and the semi norm \u\ H \ for 
u,v G V H (fl). In order not to use an unnecessary notation, we drop the subscript fl when the 
integration is over Cl and the subscript a when a = 1. 

The discrete problem associated with (1) is given by: 

Find u G V h , such that 

a h (u,v) = f{v), V v G V h (Cl). (4) 

Note that | • {h 1 (n) a nor 111 ) because if |u|^ ^(Q) = 0, then u is constant in each element. By 
the continuity at the midpoints of the edges and the zero boundary conditions, we obtain u = 0. 
Note also that / is a continuous linear form. Therefore, we can apply the Lax-Milgram theorem 
and find that there exists one and only one solution of the discrete equation (4). 

We also define the weighted L 2 norm by: 

IMIij^E) := / a(x) |u(x)| 2 dx for u G (V h + V H + L 2 )|g. (5) 
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We introduce the following notation: x < y , f >z g and u x v meaning 

x<Cy, f>cg and cv <u <Cv, respectively. 

Here C and c are positive constants independent of the variables appearing in the inequalities and 
the parameters related to meshes, spaces and, especially, the weight a(x). 



Sometimes it is more convenient to evaluate a norm of a finite element function in terms of the 
values of this function at the nodal points. By first working on a reference element and then using 
the assumption that the elements are shape regular, we obtain the following lemma: 

Lemma 1 For u € Vjg, 

H u II L a h( s ) x 5Z a CO + u 2 (M 2 ) + u 2 (M 3 )) (6) 

re r h |£ 

and 

N^ fc(E) x £ a(T) {(uiMr) - u(M 2 )) 2 (7) 

re 

+(u(M 2 ) - u(M 3 )) 2 + (u(M 3 ) - u(M!)) 2 }, 
where Mi , M 2 , are the midpoints of the edges of the triangle T as in Figure 1. 


An inverse inequality can be obtained by using only local properties. It is easy to see that for 
u£V h , 


( 8 ) 


ADDITIVE SCHWARZ SCHEMES 


We now describe the special additive Schwarz method introduced by Dryja and Widlund; see 
e.g. [11,12]. In this method, we cover Q by overlapping subregions obtained by extending each 
substructure to a larger region f^. We assume that the overlap is Si, where Si is the distance 
between the boundaries dQi and Sfl' , and we denote by S the minimum of the S,. We also assume 
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that dtl'i does not cut through any element. We make the same construction for the substructures 
that meet the boundary except that we cut off the part of fl' that is outside of fl. 

For each Qfa a P\ nonconforming finite element subdivision is inherited from the /i-mesh 
subdivision of f l. The corresponding finite element space is defined by 

V* := {u | v € V h , support of v C fi'*}, i = (9) 

The coarse space V 0 h C V h (fi) is given as the range of Ifa (or I fa) where the prolongation 
operator Ifa (or I fa) will be defined later. 


Our finite element space is represented as a sum of N + 1 subspaces 


1? 

+ 

-i- 

+ 

II 

(10) 

We introduce operators Pi : V h — ► V/ 1 , i = 0, • • • , N, by 


a h (PiW,v) = a h (w,v), VvG Vf 1 , 

(11) 

and the operator P :V h —> V h , by 


P = Pq + Pi + ■ • • + Pn. 

(12) 

In matrix notation, P 0 is given by 


P 0 = lfa(lfa T Klfa)-'lfa T K 

(13) 


where K is the global stiffness matrix associated with a/,(-, •). 


We replace the problem (4) by 

N 

Pu = g, g = Yj9i where g t - PiU. (14) 

*=o 

By construction, (4) and (14) have the same solution. We point out that gi can be computed, 
without knowledge of u, since we can find gi by solving 

a h (9i,v) = a h (u,v) = f(v), to G Vf . (15) 

The operator P is positive definite and symmetric with respect to a h (-, •). We can therefore solve 
(14) by a conjugate gradient method. In order to estimate the rate of convergence, we need to 
obtain upper and lower bounds for the spectrum of P. A lower bound is obtained by using the 
following lemma: cf. Zhang [13,14]. 

Lemma 2 Let P, be the operators defined in equation (11) and let P be given by (12). Then 

a h (P~ 1 v,v)= min ^a h (vi,Ui), Vi € Vj 1 . (16) 

Therefore, if a representation v = can be found such that 

^2a h (vi,Vi) < C$a h (v,v), to € V h , (17) 

i=0 

then 


A min(P) > Co*. 
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An upper bound on the spectrum is obtained by bounding 

a h (Pv,v ) = a h (P 0 v,v ) + a h {PiV,v ) H \-a h {PNV,v ) (18) 

from above in terms of a h (v,v). Using Schwarz’s inequality, the fact that the Pj are projections, 
and that the maximum number of regions that intersect at any point is uniformly bounded, it is 
easy to show that the spectrum of P is bounded above by 

max{#(i:p€ fl') + 1}. 
pe n 


PROPERTIES OF THE P x NONCONFORMING FINITE ELEMENT SPACE 


We first define two local equivalence maps in order to obtain some inequalities and local 
properties for our nonconforming space. Through these mappings, we can extend some results that 
are known for the piecewise linear conforming elements to our nonconforming case. 


We use a bar to denote conforming spaces. Let V 2 1^, be the conforming space of piecewise linear 
functions in Q it where the h/2-mesh is obtained by joining midpoints of the edges of elements of 

T k In,- 


We define the local equivalence map Mi : V h |^. — *■ V* as follows: 


Isomorphism 1 Given u € define u = MiU by the values of u a,t, the three sets of points (cfi. 

Figure 2.): 

i) If P is a midpoint of an edge of a triangle in T h , then 

u(P) := u(P). 


ii) If P is a vertex of an element in T h and belongs to the interior of n it and the Tj are 
the elements that have P as a vertex, then 


u(P) := mean of u^P)- 

Here u|r_, (P), is the limit value of u{x) when x € Tj approaches P. 

in) If Q is a vertex of T h \dn, , and Qi and Q r the two midpoints o/T h |an, that are next, 
neighbors ofQ, then 




Here |Q r Q| is the length of the segment Q r Q. 


Case ii) is illustrated in Figure 2., where 

O(F) = i E«lr,(F). 

° «=1 
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Case iii) is required in order to have property (21), which will be very important in our analysis. 



Lemma 3 Given u € V* 1 ^, let it € V* | ^ given by u — MiU. Then 



(19) 

NIlko.) - , 

(20) 

/ ti(s) ds= u(s ) ds. 

'dfli JdSli 

(21) 


Here | • is the standard weighted energy semi norm, for conforming functions. 

Proof. We first note that we have results similar to (6) and (7) for the conforming space V 2 , 

where now Mi, M 2 and M 3 are the vertices of a triangle in Tj. In order to prove (19), we compare 
(7) with the analogous formula for the piecewise linear conforming space. 

For instance (see Figure 2.), 

|fi(Q) - o(<MP = j^j W<3.) - «(<3.)| ! . 

The right hand side can be controlled by the energy semi norm of u restricted to the union of the 
triangles T 7 , T 8 and T 9 . 

We also prove that if we take next two neighboring vertices of 7” 2 in the interior of fij , the 
energy semi norm can be bounded locally. If a(x) does not vary a great deal, we can work with 
weighted semi norms. Using the fact that our arguments are local, it is easy to obtain the upper 
bound of (19). 

The lower bound is easy to obtain since the degrees of freedom of V h are contained in those of 

Vi 

Similar arguments can also be used to obtain (20). 

Finally, it is easy to see that (21) follows directly from iii) even if the refinement is not 
uniform. □ 
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We define another local equivalence map Mf : V h \ft. — ► V*^., by: 


Isomorphism 2 Given u £ and an edge E of dflj , define u = Atfu by the values of u at the 

three sets of points (cf. Figure 2.): 

i) Same as step i) of Isomorphism 1. 

ii) Same as step ii) of Isomorphism 1. 

Hi) If V is a vertex T h \oci, and an end point of E, and V r the midpoint ofT h \E that is 
the next neighbor of V , then 

u(V) := u(V r ). 


iv) If Q is a vertex ofT h \on x and we are not in case Hi), then 

-/vvv IQiQI /x-j \ . IQrQI •> 

U \Q) — ^ MQi) +T7TFrMQr)- 


I QlQr 


I QlQr 


Using the same ideas as in Lemma 3, we can prove: 


Lemma 4 Given u £ let u £ V 2 |q. given by u = Atfu. Then 

Mtfjcfii) x lultfi jn.) , 

IMIr-Ka) x , 

UTld 

/ iz(s) ds= u(s) ds. 

J E J E 


( 22 ) 

(23) 

(24) 


THE INTERPOLATION OPERATOR 

Let v £ V h and let Pij be the midpoint of the edge Eij common to fij and Clj. 

Definition 2 The Interpolation operator if? :V h —y V H , is given by: 

:=|^] / £)j <4,M dx = ^ 4 dx. (25) 

The second equality follows from the fact that the mean of v on each edge of an element of T h is 
equal to v(Mf), where Mi is the midpoint of the edge. It is important to note that the value of 
(if? v) (Pij) depends only on the values of v on the interface Eij. This allows us to obtain stability 
properties that are independent of the differences of a(x) across the substructure interfaces. 

Before studying the stability properties of this operator, we need two lemmas for the piecewise 
linear conforming finite element space. 

The following lemma is a Poincar 4- Friedrichs inequality. The idea of the proof can be found in 
Ciarlet (Theorem 6.1) [9] and in Necas (Chapter 2.7.2) [15]. 
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(26) 


Lemma 5 Let T be a subset of dQi, such that T and d£l{ have measures of order H. Then, 

II^IIl 2 ^) — Vue H 

As a consequence, if J r u(x)dx = 0, we have the Poincare inequality 

||w||L 2 (n< ) r< H lul/rjtni)- 


(27) 


The next lemma is a Poincare-Friedrichs inequality for nonconforming Pi elements. It is obtained 
by using Lemmas 3, 4 and 5. 


Lemma 6 Let u G Hl h (Q,i), where f2j is a triangular substructure of diam.eterO(H). Let T be dtl, 
(or an edge of d£li). Then, 

Nil Hai) r<W\u\l lw + (j T u(x)dx)\ VueHtia,). (28) 

As a consequence, if f r u(x) dx = 0, we have the Poincare inequality 


Nlr£ h (fti) ^ H MtfVfii)- 


(29) 


The next lemma gives an example of an operator that is L^— and —stable. 

Lemma 7 Let u G tf*(Qi), where is a triangular substructure of diameter of 0(H). Define a 
linear function Ujj in flj by 



u H (Pij) [ H x ) dx > j = l,2,3, 

\btij\ JEij 

(30) 

where the Eij 

are the edges ofQ.i, and Pij is the midpoint of Eij . Then , 



|WH(-P»i)| 2 d -jp\M\LH(li) + 

(31) 

and 

luffing) ^ 1 “lifKa). 

(32) 

II^h - ^IU 2 (n,) d # N^n,)- 

(33) 


Proof. Consider initially a region S7 with a diameter of 1. Using that \Eij\ — 0(1), the 
Cauchy-Schwarz inequality and a trace theorem, we have 

|0 H (Py)| ! x I / 0(z) ^ ll“lli=(E„) 

JEij 

- - ll fi H Wo) + 1*1 HHEij)- 

We obtain (31) by returning to a region of diameter H. 
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Note that for any constant c 


(34) 


|w/f(-f»l) — Wh(-F» 2)| 2 + |Wi/(-F»2) — Wtf(-P»3)| 2 + \un(Pi3) ~ itff(-Ftl)| 2 

d 11^ — c llH 1 (O i )- 

By choosing c = u(Pn) and V = £?«., we can apply Lemma 5 and obtain the H 1 -stability (32). 

We now prove the Instability. Since u — un has mean zero on di li, we can apply the Poincar6 
inequality (27) and obtain 

ll« — d H\u- unlimsii)- 

Using the first part of this lemma, we obtain the L 2 -stability (33). □ 

The next lemma shows that the interpolation operator Iff , defined by (25), is locally L\— and 
H\— stable. 

Lemma 8 Let u € V h (Cl). Then Uh — Iffu satisfies the following properties 

Wh\hi ff (0j) d l u l«2 ifc (a-) * ( 36 ) 

CLTtd 

||«h — u IU|(n<) d H |u| H i h ( n .) , i = ,N. (37) 

Proof. Let uh = Iffu and let u € H 1 ( fi*) be given by u = A if n u and let uh(Ph) be given by 
(30). Using the properties (24) and (25), we have 


u H (Pi i) = u H (Pn). (38) 

Therefore, by (38), (31) and Lemma 4, we have 

Wh(Ph )\ 2 = |«/f(P»i)| 2 d jp\\u\\h{{ii) + ( 39 ) 

d + Mnqn,-)- 

We also obtain the same estimate for |uh (P^)! and |utf(P i3 )|. 

The rest of the proof is similar to that of Lemma 7. We now use the Poincare inequality for 
nonconforming elements. □ 


THE PROLONGATION OPERATOR 

In this section, we introduce several prolongation operators and establish that they are stable. 
The range of each of these operators will serve as a coarse space in our algorithms. 
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Definition 3 The Prolongation Operator 1% : V H -* V h , is given by: 

i) For all nodal points P ofT h that belongs to an edge Eij common to Q* and Clj, let 

{I H u h){P ) := where Pij is the midpoint of the edge Eij. 

ii) Given IjjUh at the nodal points of T = Ujdfli from, i), let be the 

P\-nonconforming harmonic extension inside each ■ 

It is easy to check that = I^uu £ V h {fit). A disadvantage of step it) is that we have to solve 
exactly a local Dirichlet problem for each substructure in order to obtain the harmonic extension. 
Other extensions can be used, which we call approximate harmonic extensions. They are given by 
simple explicit formulas and have the same L\ and H\ h stability properties as the harmonic one. 



Our first construction is a natural generalization of the partition of unity introduced by Dryja 
and Widlund in [11]; this partition of unity will provide the basis functions of our approximate 
extensions. Let Pj, j = 1, 2, 3, be the midpoints of the edges of 0*, and let Vj be the vertex of 
that is opposite to Pj. Let C be the barycenter of the triangle fli, i.e. the intersection of the line 
segment connecting Vj to Pj. 

Extension 1 The construction of an approximate harmonic extension is defined by the following 
steps (see Figure 3.): 

i) Let 

u(C) := ^ { u H (Pi ) + ii# (P 2 ) + uhCPs)}. 

ii ) For a point R that belongs to a line segment that connects C to a vertex Vj , let 

u(R) := u{C ). 

Hi) For a point Q that belongs to a line segment connecting C to Pj, define u(Q ) by 
linear interpolation between the values u(C) and Uu{Pj). i- e by 

u(Q) := A (Q)u(C) + (1 - A (Q))u H (Pj). 


Here A (Q) =dist,ance(Q, Pj)/ dist,ance(C, Pj) . 
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iv) For a point S that belongs to the line segment connecting the previous point Q to a 
vertex V*, with k ^ j, let 

u(S ) := Ci(<?). 

v) Finally, let I^Ufj = I^u, where Ih is the interpolation operator into the space V h that 
preserves the values of a function at the midpoints of the edges of the elements. 

Note that the function u just constructed is continuous except at the vertices V, of fl*. The step 
i) can be viewed as emulating the mean value theorem for harmonic functions. However, near the 
vertices, u is a bad approximation of the harmonic extension. We know that the local behavior of 
the harmonic extension near a vertex Vj depends primarily on the boundary values in the vicinity 
of Vj. For instance, if ujj{P\) = 0, uh{Pz ) = 0, and uh{Pi) = 1 , we should obtain Uh — 0 near Vi ; in 
addition, by using symmetry arguments, we should have u h ~l/2 for points near Vi that lie on the 
bisector that passes through Vi . With this in mind, we now construct an alternative approximate 
harmonic extension. 

We change notation in order to be able to use Figure 3. Now let C be the point where the three 
bisectors intersect. 

Extension 2 The construction of the approximate harmonic extension is defined by (see Figure 

3.): 


i) Same as Step i) of Extension 1. 

H) Define u(Vj ) — | u(Pi) ■ Dor a point R that belongs to a line segment connecting 
C to Vj, define u(R ) by linear interpolation between the values u(C ) and tt(Vj). 

Hi) Same as Step Hi) of Extension 1. 

iv) For a point S that belongs to a line segment connecting the previous point Q to 
Vic, k j, u(S) is defined by linear interpolation between the values u(Q ) at Q and 
f(QJ,k) at Vk. Here, 


f(Q,j,k) = MQ)u(V k ) + (1 - A (Q))u(Pj). 


v) Same as Step v) of Extension 1. 


A disadvantage of this extension is that we cannot just work in a reference triangle, since the 
angles are not preserved under a linear transformation. This is similar to the fact that under a 
linear transformation a harmonic function does not necessarily remain harmonic. We can construct 
other approximate harmonic extensions which combine the properties of the two extensions, given 
so far, and working, for instance, with the barycenter C as in Extension 2 and replacing the weight 
1/2 in Step ii). 

The next lemma shows that the extensions given above have quasi-optimal energy stability. 
Using ideas of Dryja and Widlund[ll], we prove the following lemma. 
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Lemma 9 Let uh € V H (Q,). Then 

\ j h u h\hI h (fii) r< (1 + log(H/h))~- \u h \ h i (40) 

and 

\\Ih u h ~ uhWlkq,) d H |ittf|tfi w (n,)- ( 41 ) 

Proof. Let 0 3 h € V h \ n 0 j = 1,2,3, be the approximate harmonic extensions constructed from the 
boundary values = 1 at the /i-mesh nodes on the edge Eij, and 0^ = 0 at the other boundary 
nodes of dQ,i. It easy to see that the 6 3 h form a basis of all approximate harmonic extensions that 
take constant values on the edges of the substructure. It is easy to show that if a point x belongs to 
the interior of an element of fij, then |V #^(x)| is bounded by C/r, where r is the minimum distance 
from x to any vertex of 0*. Note that any element that touches a vertex of provides an order one 
contribution to the energy semi norm. To estimate the contribution to the energy semi norm from 
the rest of the substructure, we introduce polar coordinate systems centered at the vertices of 0,%. 
Then, 

|0£|^ (n .) - 1 + / J h r ~ 2 r dr d( P ^ 1 + log(H/h). (42) 

Since the partition of unity forms a basis, it is easy to see that 

^ ( 43 ) 

(1 + log(H/h )) {|uh(Ti)| 2 + \uh{P 2 )\‘ 2 + |un(P3)| 2 } 
and using ideas similar to that of Lemma 7, we have 

^ (! + log(H/h)) flu# (Pi) - u H (P 2 )\ 2 + 

|utf (P 2 ) — Uh(P3)\ 2 + ^//(Pj) — Uh(P\)\ 2 } 

x (1 + log(H/h))\uH\ 2 H i^ n .y 


By construction, it is easy to see that 

\(I h H u H )(x)\ < ,max 3 |iijf(P*)|- 

Therefore 

II I h „UH - «Hlll !( n,) iE"’ M«)l’ . 

i 

and by using (39) and (29), we obtain (41). 

Since a(x) varies little in each these arguments are also valid for the weighted norms and we 
obtain (40). □ 
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Using Lemmas 6 and 9 and the triangular inequality, we have: 

Theorem 1 Let u € V h (Q). Then 

Wh 1 ^ u ~ w IUi(n<) ^ H Mtf] h (tti) (44) 

and 

l^H^h u \h\ ^ (1 + log{H/h)f \u\ H i (45) 

Remark 1 It is easy to see that we do not need to use the fact that ujj 6 V//(51); we only need to 
calculate values Vfj(Pij) by formula (25) at the midpoint P^j of the edge Eij. The next step is to 
provide the constant value Vn(Pij) to all nodes of the interface and perform, an approximate 
harmonic extension. 

Remark 2 The extensions also can be constructed for nontriangular substructures. In a, first, step, 
we construct a partition of unity in fi*. This can be done by using ideas similar to those of the 
triangular case. By using the same technique as in the proof of Lemma 9, we can show that 

\Ih u h\hI h (cii) ^ (46) 

(1 + log(H/h)) £>(0,) MPy) - MP jU -„)| 2 

!-■ ■ 

where and P\{j- 1 ) are neighboring midpoints of edges of d£li and N* is the number of edges of 

d£li. We obtain (44) by noting that each term of the sum, is bounded by \u\ 2 H i 

a , h * 1 ' 


THE NEUMANN-NEUMANN BASIS 

In this section, we consider a Neumann-Neumann coarse space. This is the P\ nonconforming 
version of a coarse space studied in Dryja and Widlund [16], and Mandel and Brezina [17]. 

However, here we use an approximate harmonic extension inside the substructures. We note that 
the coarse spaces considered by these authors differ only in how certain weights are chosen. Mandel 
and Brezina use weights that are convex combinations of the coefficient a(x), while Dryja and 
Widlund use a?(x). Here we show that any convex combination of a /3 (x), for (3 > 1/2, leads to 
stability. We point out that the choice (3 — 1/2 can be viewed as a L 2 -average, while (3 = 1 is an 
average in the L 1 sense. 

We call the coarse space of the previous section, face based. There are some differences between 
Neumann-Neumann and face based coarse spaces. A Neumann-Neumann coarse space has one 
degree of freedom per substructure, while a face based uses one degree of freedom per edge. A 
Neumann-Neumann basis function associated with the substructure fl*, has support in f2j and its 
neighboring substructures, while a face based function basis, associated with an edge of a 
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substructure, has support in just two substructures. The face based coarse space appears to be 
more stable since all the estimates, related to the jumps of the coefficients, are tight. In the lemmas 
that we have proved for the face based methods, all the stability results were derived in individual 
substructures, while in the Neumann-Neumann case, we need to work in an extended subdomain. 



i) For each substructure calculate the mean value on dQ, if i.e . 

m ‘ u:= iakj L u{s)ds - 


v h , as follows: 


Here is the length size of dfli. 

ii) For all nodal points P ofT h that belong to the edge E iti , let 
(I NN U)(P) = (/?!»)(/>„), where 




:= miU 


Here Pij is the midpoint of the edge Eij . 

Hi) Perform an approximate harmonic extension to define InnU inside the substructures. 


Note that we can also calculate m^tt by: 

(47) 

Therefore, there exists a linear transformation : Vh —* Vh, such that Iffu = I^lffu. The next 
lemma establishes stability properties for 1%. 
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(48) 


Lemma 10 Let uh € V H (Q) and (3 > 1/2. Then 

\Ih u h\hI w (n ( ) < C{P) \ u n\Hl H (^ a, ) > 

and 

II Ih u h ~ whIUko,) < C{0) H\u h \ h \ ff (n"')- 
Here the extended domain fi® xt is the union of f 2* and the substructures that share 


Proof. Let us first prove the L\ stability. Note that (see Figure 4.) 


K(Ptf) - (/f ««)(Py) I 2 = MPij) 


a,P(£li) rrii + a^(Qj) m # 

a^(a) + ‘ 


By using (47) and simple calculations, this quantity is equal to 


1 

|a0(f2 i ) + a' J (O,)| 2 * 

!«?(«.) {j^|(K»(Py) - u„(P ik )) +j^|jW(P i j) - u„(P„))}+ 


Using the shape regularity of the subdomains, it is easy to see that 

am iu„(p y ) - asmiPti) i 2 d 

and using the fact that f3 > 1/2, we can bound this quantity by 




— C(f3) \uh\ 2 hi ^(fi.un,)- 

We obtain (49) by adding all the contributions (50) to the L 2 a (Q.i) norm. 

We prove (48) by using the triangular inequality, an inverse inequality, and (49). 
Theorem 2 Let u € V h ((l) and (3 > 1/2. Then 

HLvwu - u||^2 { n.) < C((3 ) H |u| H i h (n?*’) > 

and 

\Innu\hi ^(n,.) < C{(3)( 1 + log{H/hyp |u| H< i fc (nr")- 


(49) 

edge with f2j. 


(50) 


□ 

(51) 

(52) 
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Proof. Using Lemmas 9, 10 and 8, we have 

\Innu\hi h (fii) ^ (1 + log(H/h))* \I^I^u\ H i w (n,) < 

C(t 3) (1 + log(H/h))i |/f u|„, < 

C(0) (1 + log(H /h))- l«l»; 

The Instability is obtained by 

\\I NN u - w||i2(n,) < \\Innu ~ lH I h u \\Li(n i )+ 

11 Ih i h u - IhUWma,) + \\IhU- u|| L ;2 (nj) , 
and by using Lemmas 9, 10 and 8. □ 

Remark 3 We can also prove Theorem 2 for the case of nontriangular substructures; c.f. Bemmks 
1 and 2. 


THE THREE DIMENSIONAL CASE 


We show in this section that the methods developed before can be extended to three dimensions. 

For simplicity, we assume that is a polyhedral region of diameter 1 in three dimensional space. 
As before, we introduce a nonoverlapping partition composed of tetrahedra fli of diameter of order 
H. This defines a coarse space and a triangulation T H . We further subdivide the substructures into 
tetrahedra which results in a triangulation T h and define the nonconforming P\ finite element 
spaces V h and V H as in Definition 1. Here, the continuity is enforced at the barycenter of the faces 
of the triangulations. 

The local equivalence maps are given by the following procedure. In each tetrahedral element of 
T h (cf. Figure 5.), we connect its centroid to the four vertices and to the barycenters of the four 
faces. We also connect each barycenter to the three vertices. In other words, we subdiyide each 
tetrahedral element into twelve subtetrahedra. We denote this new triangulation by T h . The 
vertices of T h are the vertices, barycenters, and centroids of the elements of T h . 

Let be the conforming space of piecewise linear functions of the triangulation T h \ci r 

We define the local equivalence map Adj : V h — ► V h |jy , as follows: 

Isomorphism 3 Given u € V h \fi., define u = MiU by the values of u a.t. the following sets of 
points: 
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i) If P is a vertex of an element ofT h and belongs to the interior offy, and the K 3 are 
the elements in T h \^. that have P as a vertex, then 

u{P) := mean of u|j*\(P). 

Here u\ Ki {P) is the limit value of u(x) when x £ Kj approaches P. 

ii) If P is a barycenter of a triangle in T h U, then 

u(P) := u(P). 


Hi) If P is a vertex of a triangle in T h \ d n . and Tj, j = 1, • • • , N P , are the triangles of 
P h \dCi, that have P as a vertex, then 


u 


(P) 

fc=l I u j=1 lj I 


Here C, and |Tj[ are the barycenter and the area of the triangle Tj, respectively. 

It is easy to check that the Lemma 3 holds, if we replace V h ^ 2 \^i, by V^|n . 

We define another local equivalence map Mf : , by: 

Isomorphism 4 Given u £ and a face F of dQ. t , define u = Adfu by the values of u at the 

following sets of points: 

i) Same as step i) of Isomorphism 3. 

ii) Same as step ii) of Isomorphism. 3. 

Hi) Let P be a vertex of a triangle in T h \aci { that belongs to dF, and let. Tj, 
j = Np, be the triangles of T h \ F that, have P as a vert, ex. Then 


NP 


I T k 


*(P) := £ - 

fc =i I U i= p ! Tj\ 


uiCJ. 


iv) Let P be a vertex of a triangle in T h \ d Q i that does not. belong to dF, and let. Tj, 
j = 1, •• • , N P , be the triangles ofT h \ F that have P as a vertex. Then 


N P 

u(P) := Z 


|T*| 


' i \ N P nn 
k = 1 u j=l 1 j 


u(Ci). 


It is easy to check that Lemma 4 holds, if we replace V h/2 \ q. by and let the faces play the 

role previously played by the edges. 

Let v £ V h and let Cij be the barycenter of the face F {j common to Cl t and f \j. 


560 


Definition 5 The interpolation operator I ff : V h — * V H , is given by: 

(/f»)(c«) :-j^j 

where \Fij | is the area of the face ho- 
using the same ideas as in two dimensions, we can prove lemmas analogous to Lemmas 5-8. 

The prolongation operator I'fr : V H -* V h , is defined as in the two dimensional case. In a first 
step, we define (/^u H )(P) := u H {Cij) for all barycenters P of triangles in T h \ Fij . Finally, we 
perform a Pi -nonconforming harmonic or approximate harmonic extension. 

We describe the three dimensional version of Extension 1. This is a generalization of the 
partition of unity introduced by Dryja, Smith, and Widlund [14]. Let Cj, j — 1, • • • , 4, be the 
barycenters of the faces Fj of d£l % , and let Vj be the vertex of that is opposite to Cj. Let C the 
centroid of SI*, i.e. the intersection of the line segments connecting the Vj to the Cj. Let E jk , 
k = 1, 2, 3 , be the edges of dFj. 

Extension 3 The construction of an approximate harmonic extension I%Ufj is defined by the 
following steps (see Figure 5.): 

i) Let 

O(C) := i X>«(CA 

4 i= 1 

ii) For a point Q that belongs to a line segment, connecting C to Cj, define u(Q) by linear 
interpolation between the values u(C) and ujj{Cj), i.e. by 

u{Q) := A (Q)u{C) + (1 - A(Q))uh(C , j ). 


Here A(Q) =distance(Q,Cj)/distance(C,Cj). 

Hi) For a point S that belongs to any of the three triangles defined by the previous Q, and 
the edges Ej k , k = 1, • • • , 3, let 

u(S) := u(Q). 

iv) Finally, let I^u H = hu, where I h is the interpolation operator into the space V h that 
preserves the values of a function at the bary center of the faces of elements in T . 


We can also construct an approximate harmonic extension similar to that of Extension 2. This 
gives a better approximate harmonic extension near the edges. 

The prolongation operator Ifj in three dimensions has the same stability properties as in the two 
dimensional case, i.e. Lemma 9 still holds. 


561 


The idea of the proof is the following. Consider the case where is given by uh(Ph) = 1 

and u H {P i2 ) = un{Piz) = 0. This gives the partition of the unity introduced by Dryja, Smith, and 
Widlund [4]. The energy semi norm of uh is of order H. 

Let 0 j, 1 = We note that is bounded by C/r, where r is the distance to the 

nearest edge of f The contribution to the energy semi norm from the union of the elements with 
at least one vertex on the edge of the substructure can be bounded by CH , since the extension is 
given by a convex combination of the boundary values. To estimate the contribution to the energy 
from the rest of the substructure, we introduce cylindrical coordinates using the appropriate 
substructure edge as the 2 -axis. Integrating |V0J l 1 (x)| 2 over this region, we find that it is bounded 
by C (1 -I- log(Hfh)) H. 

To prove Lemma 9 for a general uh, we use the same ideas as for two dimensions. Similarly, we 
can extend the results to nontriangular substructures and to the Neumann-Neumann case. 



MAIN RESULT 

In this section, we consider the Schwarz method introduced in the previous sections and prove 
the following result. 

Theorem 3 The operator P of the additive Schwarz algorithm., defined by the spaces V 0 h <ii id v t h , 
satisfies: 

K (p)i(i + ; 09 (f))(i + f)- 
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Here k(P ) is the condition number of P. Therefore, if we use a generous overlapping , then 

k{P) i< l + log(H). 

Proof. The proof of this theorem is essentially the same as in the case of a conforming space; see 
Dryja and Widlund [12]. 

As we have seen before, the upper bound is very easy to obtain. The lower bound is obtained by 
using Lemma 2. We partition the finite element function u 6 Vh as follows. We first choose 
uo = Mu or Innu, i.e. apply a face based or Neumann-Neumann interpolation operator. Let 
w = u — uo- The other terms in the representation of u are defined by u, — i = 1, • • • , N. 

Here J/, is the linear interpolation operator into the space V h that preserves the values at the 
midpoints of the edges of the elements and {#i} is a partition of unity with 9, € C$)°(f2i) and 
£0i(*) = i. 

For a relatively generous overlap of the subdomains, these functions can be chosen so that V 0, is 
bounded by C/H. By using the linearity of Ih, we can show that we have a correct partition of u. 
In order to estimate the semi norm of Ui, we work on one element A” at a time. We obtain 

JK) — 2 \^i w \'hI h (K) + 2 |Jfc((0* — Oi) w )\Hf h {K) 

Here §i is the average value of 9, over K. It is easy to see, by using the inverse inequality (8), that 

\h({9i - Oi)w) \]ji ■< h~ 2 ||I/ l ((0i - 9i)w)\\\2^ K y 

We can now use the fact that on K, 6, differs from its average by at most Ch/H. After 
s umming over all elements of we arrive at the inequality 

Nifi jn') ^ klHi h (n;.) + H ~ 2 \\ w \\hw 

We sum over all i and use that each point in is covered only a fixed number of times and 
obtain a uniform bound on Cq. We conclude the proof by estimating the two terms of 

Mtfi h (n) + H 2 IMIl;j(fi) 

by |u|^i The bounds follow by using the stability results of Theorem 1 or 2. 

For the case of small overlap, the proof is similar to that of the case of piecewise linear 
conforming space considered in Dryja and Widlund [12]. 0 
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